Study of B° -> D°ji + ii Decay in Perturbative QCD Approach 
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Abstract 

Within the perturbative QCD approach and ignoring the contributions of long distance and subleading penguin loops, 
we investigate B° — > D°ii + ii~ decay in the large recoiling kinematic region in the Standard Model. At the tree level, 
B° decays to D° by exchanging a W boson accompanied by a virtual photon emission from the valence quarks of 
B° and D° meson, then the virtual photon decays to the lepton pair. Numerically, we find that the branching ratio 
decreases rapidly as the q 2 increases, and the branching ratio of B° — > D° is (9.7+^) x 1(T 6 in the region 

q 2 e [1 , 5] GeV 2 . The order of the branching ratio shows a possibility to study this interesting channel in the current B 
factories and the Large Hadron Collider. The precise experimental data will help us to test the factorization approach 
and the QCD theory, in general. 
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Over the past few years when studying the semileptonic decays of B meson, people always pay much attention 
on exclusive processes B — > (K,K* , 7t,p)£ + £~ and inclusive processes B —> X s ^t + t~ as well as similar decay modes, 
which are induced by the flavor changing neutral current b — > si + £~ or b — > d£ + £~ . In these processes, the leptons are 
always generated from either a photon or a Z boson with loop diagrams, so that these decay processes are considered 
as good choices of testing the Standard Model (SM) and probing possible new physics signals. Recent review in 
detail is referred to Refs. fU|2][3|. In fact when we study the decays B — > (K,K* ,7t,p)£ + £~ , the weak annihilation 
contributions are usually ignored since they are regarded to be suppressed by <^(Agc£>/ m fl) B- Therefore, we think 
that it is of urgent interest to explore the pure annihilation type semileptonic B meson decays, in which ^(AgcoAwfl) 
effects are the main contribution. Still due to suppression of ^(AgcD/mg), most of these decays have small branching 
ratios, and cannot be observed in the current BaBar and Belle experiments. However, for some special decays, such 
as B° — > its branching ratio can be enhanced by large Wilson coefficients. In this work, we consider the 

observables of the decay B° —> D jj, + jj,~ theoretically. Compared with the mass of B meson, both the masses of muon 
and electron are very small, so the analysis of B° — > D°e + e~ is almost the same as 

In the SM for B° £>W~ the muon pair can be generated from either a photon or a Z boson, however, the 
latter case will be highly suppressed because of the weak coupling and the large Z mass. Therefore, we only consider 
the process where the lepton pair is generated from a virtual photon. In the full theory, there are three possible 
contributions to this decay, and we draw the Feynman diagrams in Fig. [TJ In the first case, shown in diagram 1(a), B° 
decays to D° + J / iff by exchanging a W boson and generating cc pair from the vacuum, in which the J / y/ decays to 
lepton pair, which is so-called the resonant contribution. Because the mode B° — > D° +//!//■(—> £ + £~) has not been 
observed yet, we will exclude this part of contribution, i.e. Fig. 1(a), by carrying out our investigation in a certain 
kinematics region, q 2 <C m ^j/y The virtual photon can also be generated by the penguin operator or Ofjy, which is 
shown in diagram 1(b), with the Wilson coefficient Cj. Since this operator is from the loop suppressed flavor changing 
neutral current, the value of C7 is much smaller than those of the coefficients C1.2 of tree operators, and thus only 
marginally affecting our numerical estimates. Therefore, the contribution of diagram 1(b) has been neglected safely 
in this work. In diagram 1(c), the B meson decays to a D meson by exchanging a W boson, where the photon can be 
emitted from either of the five crosses in diagram. When a photon is emitted from the W boson, the diagram will be 
highly suppressed by the two W propagators and because of the large W mass. Therefore, we ignore this contribution 
in our calculation, too. Since this process happens at the scale ^(mj), the highly off-shelled W boson can be integrated 
out and the effective theory could be used directly, as shown in Fig. 2. 
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Figure 1 : The possible diagrams for B° -> D°£+£-, where the crosses stand for a virtual photon. 
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To make predictions clear, one requires the knowledge of the matrix element (Dy* \B), where the virtual photon y* 
decays to a lepton pair. Although the calculation of this matrix is not trivial, it has been explored in many approaches, 
such as the heavy quark effective theory the heavy light chiral perturbation theory [6|, the QCD factorization 
approach Q and the perturbative QCD (pQCD) approach (8). Based on kj factorization, the pQCD approach ll9l [T0l is 
one of the theoretical instruments for handling such exclusive decay modes. The concept of pQCD is the factorization 
between soft and hard dynamics. In this approach, the quark transverse momentum kj is kept in order to eliminate the 
end-point singularity. Because of inclusion of transverse momenta, double logarithms from the overlap of two types 
of infrared divergences, soft and collinear, are generated in radiative corrections. The resummation of these double 
logarithms leads to a Sudakov factor, which suppresses the long-distance contribution. Though there still exist few 
controversies |[TT1|12I on its feasibility, the predictions based on the pQCD can accommodate experimental data well, 
for example, see Ref. 021. In this work, we will put the controversies aside and adopt this approach to our analysis. 

In the SM, the effective Hamiltonian related to decay B° — > D°£ + £~ is given lfl4l as: 



j% s = ^Jf cb v: d [C, (m) +C 2 ( M )0 2 (m)] , 



(1) 



where Gf is the Fermi constant and V c t,V* d are the corresponding CKM matrix elements. Oi and O2 are local operators, 
which are defined as: 



0\ — {c a bp)v-A{dpu a )v-A , 

2 = {c a ba)v-A(dpUl})v-A ■ 



(2) 



Here a, j3 are the color indices, {q\q2)v-A = <?i (1 — J 5 )q2, and C\ and C2 are corresponding Wilson coefficients, 
whose scale evolves from mw to the factorization scale t . With the Hamiltonian in Eq. (fl}, we draw the diagram in 

Fig.m 

Now, we turn to discuss the decay B° — > D°jx + n~ in certain kinematic regions like V stalt < q 2 < V ea d, where q is 
the momentum of the £ + £~ pair, y sta rt and V en( j are the boundaries of the region. To guarantee our calculation reliable, 
we should choose the region where D meson recoils fast and it can be treated on or nearly on the light cone. In the rest 
frame of B meson, the momenta of B and D mesons are defined in the light-cone coordinate as 



m B 



= ^(i,i,3 ± ) , pd = ^|(t? + vV-i^- vV-1,60 



(3) 




(a) 




(6) 



Figure 2: Diagram for B° — > D£ + £ in the effective theory. The black boxes represent the effective vertex. 
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with 

m- B + m- D - V end 



1 mn-V e „ c \ ml+mn-V s t 



For the light quarks in B and D mesons, we define their momenta as 



<?7 < 5 - u — . (4) 



l !(\ nlB 7 \ 1 / r l + \/ r l 2 ^ 1 a 7 \ 

x 2 m D ,0,k 2 j_) , (5) 

where k± stands for the transverse momentum. 

For the decay B° — > D°£ + £~ the amplitude will be factorized conventionally to a hadronic part and an electromag- 
netic part. To make our expressions simple, we parameterize the hadronic matrix element with two contracted weak 
vertices and one QED vertex as 

= (D \Ci(n)Oi(n) ^ \B ) =fi{q )p%+f 2 {q )Pu , (6) 

where f\ (q 2 ) and fiiq 2 ) are form factors, and their expressions are given by 

Mq 2 ) = fiA(q 2 )+fi.2(q 2 )+fiAq 2 )+fiA(q 2 ), 

fl(q 2 ) = f2,l(q 2 )+f2,2(q 2 )+f23(q 2 )+f2A(q 2 ), (7) 

in which the second subscripts of f l} correspond to the numbers of the crosses in Fig. |2] Within the perturbative QCD 
approach, in the large recoiling region, fi j could be calculated at the leading order up to the leading power of mo / ms. 
The detailed expressions are given in Appendix lAl Unlike the form factors of the charged current process B° — > D , 
fi and f2 are complex numbers, which are caused by the annihilation mechanism. Numerical results in the region 
q 2 £ [lGeV 2 ,5GeV 2 ] show that both the real and imaginary parts of f\ are much larger than those of fi. 
With the functions defined above, the amplitude can be expressed as 

Jt= C jf cb V* d T»[lY ll l] =Mq 2 )[l fal}+f 2 (q 2 )[l Pol], (8) 



and 



with 



S\\ = Tr[{h + m i) Mfa-mi) Pb], 
Si2=Tr[(fii+mt) j) B {fa-mi) p D ], 
S2\=Tr[{p l +m,) p D (p 2 -mi) Pb], 

S22 = Tr[( Pl +m l ) Mh-m) M (10) 

In the above functions, p\ and p 2 are the momenta of the l~ and / + leptons respectively, and m/ is the lepton mass. In 
the center of mass frame for the lepton pair, we define p\ and p' 2 as corresponding momenta of p\ and p2, 

P\ = (v^/2,psin0cos0,/?sin0sin0,/7cos0), 

p' 2 = {\fq 2 ~/2,— psin0cos0,— psin0sin0,— pcosO), (11) 
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where p is the magnitude of 3-component momentum and p 2 = q 2 /4 — mf, 9 [0] is the inclination [azimuth] coordinate 
of l~ . After the Lorentz transformation, one can get the expressions for p\ and p2 as follows. 

Pi — (yv?/2 — 7/3/7 cos 0,/7 sin 0cos0,/?sin 9 sin0, — 7j3 \fq 2 /2 + ypcosQ), 

P2 = (7\/?/2 + 7j3pcos0, — psin0cos0, — psinO sin</>, — yfi yfq 2 /2 — 7/?cos0), (12) 

where f5 = "]^_Jl~ D ^ and 7 = (1 — /3 2 )~ 1//2 . As a consequence, the expressions for Sy with i,j = 1,2 are given as 

5n = m 2 B (4m 2 cos 2 9 + q 2 sm 2 9)[-l+Y 2 (l+l5 2 )] 1 



Sn 

521 
^22 



OT B m D (4m 2 cos 2 9 + q 2 sin 2 9) 



-tj + 777 2 (1 +/3 2 )+2j37Vi7 2 -l 



^(4m 2 cos 2 0+^ 2 sin 2 0)|-l + 7 2 -1 +2r/ 2 + j8 2 (2r/ 2 - l) +4/3r/ vV - 1 }. (13) 



The most important inputs of the calculation are hadron distribution amplitudes, named <ps and <pp, which contain 
the nonperturbative effects in the mesons under the scale Aqcd- Under the factorization frame, they are universal 
quantities and can be constrained from well measured decay channels. For the B meson distribution amplitude, we 
adopt the model [9 1 : 



$B(x,b) =Nbx 2 (1 -x) 2 exp 



2a2 



_1 (xM B y co 2 b 
2 \ cob J 2 



(14) 



with the shape parameter cob = 0.40 ±0.05 GeV, which has been tested in many channels such as B — !► nn,Kn 110]. 
The normalization constant is related to the decay constant fg = 190 MeV J9) by the normalization condition in 
Eq. ( [Tol l. As for D meson, the distribution amplitude, determined in Ref. (TBI by fitting, is 



0D 



2V6 



f D 6x(l -x) [1 +C D (l -2x)]exp 



orb 2 



where Cd = 0.5, ft) = 0.1. Both distribution amplitudes are normalized as: 



Jo 



dx<p M (x) 



fi 



M 



2V2N^ 

One can obtain the differential decay width by 

dT _ VI 
dq 2 d cos 9 dtp ~ m2An A m\ \ 



M = B,D. 



4mj 



(15) 



(16) 



(17) 



where X = (m B + nr D — q ) — \nr B nr D . Integrating over the angle variables, we would obtain the q -dependance of 

I 2 _ 4m 2 

the decay width as well as the branching ratio. In Eq. ( fT7| ), the factor ■ ' - — — 



ensures that the branching ratio 
at q l = 4mf vanishes, however, the q 1 appearing in the denominator of the photon propagator generates a pole-like 
structure at the small q 2 region. Since it is very difficult for the detector to observe leptons with such a low energy, 
we simply subtract the region with very small q 2 value. In addition, in order to avoid the pollution from long distance 
contributions shown in Fig 1(a), we set the maximum value of q 2 as 5 GeV 2 . 

In Fig. [3] we present the behavior of the branching ratio of this decay mode with 1 GeV 2 < q 2 < 5 GeV 2 . From 
the figure, one can see that the value of the branching ratio decreases rapidly as the q 2 increases: at q 2 = 1 GeV 2 the 
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Figure 3: The dependence of the branching ratio of B° — > D°jj. + jj. with q 2 , and q 2 € [1 , 5]GeV 2 . 

value is 3.2 x 10~ 5 , and it decreases to 2.8 x 10~ 8 at q 2 = 5 GeV 2 . By integrating the branching ratio over q 2 in the 
region [1,5] GeV 2 , we obtain: 



where the errors are mainly from Aqcd- The errors from the decay constant are not listed directly, which are propor- 
tional to the square of the decay constants. We here do not discuss the uncertainties taken by CKM elements, simply 
because they can be measured well in other decay channels. Since there only vector currents appear in the calculation, 
there is no forward-backward asymmetry in this decay mode at the tree level, so any apparent deviation from zero 
would be the signal from new physics. The order of magnitude for branching ratio shows a possibility to study this 
channel in present Belle, BaBar and LHC-b as well as future Super-B factories. The precise experimental data will 
help us to test the factorization approach, and the QCD theory itself in general. We are pretty sure that future studies 
on the decays will come soon from several other theoretical approaches, and the numerical estimates will be further 
refined. 

Finally, let us summarize our work. Within the pQCD approach, we studied the exclusive rare decay of B° — > 
D°n + [i~, which is pure annihilation type decay. Explicitly, we have found that the branching ratio is (9.7^3 2) x 10~ 6 
and the forward-backward asymmetry is zero at the tree level. It is clear that such an order of magnitude for branching 
ratio could be well measured at the ongoing B factories and Large Hadron Collider as well as future Super-Z? factories. 
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A Appendix A: Relevant Functions 

The definitions of fu used in the text are presented in this appendix. These functions can be calculated directly within 
the perturbative QCD approach: 

2 f l f A QCD 

/u(<7") = 4e b 7ia em m D f D / dx\ / diifeia 2 (fi)exp[-5 B (fi)]fe(xi) 7 -2 — — — H^D x b x ) 



o JO " (r\--\)q 2 n 



x ( 2m D (4?7 4 - 5?7 2 - 3T] ^?7 2 - 1 +477 Y?] 2 - 1 + 1) - m B (x! - 2)(2?7 3 + 27] V?] 2 - 1 - ^/t] 2 - 1 - 27]) ) , 

AlW ) = 4e b 7l(X em m B f D / rfxi / rf&iZ7 1 a 2 (fl)exp[-5B(fi)]0 B (jci) 7 -2 — 7y^- /:/ o(\ADl£'i ) 

JO {r\--\)q A K 



m B (x 1 - 2)(7 7 2 + 77 ^77 2 -l-l) + 2 mD (-277 3 -277y77 2 -l + ^7 7 2 - I+277) ), 

1 , /' Aqcd ..x r „ /. sn j. / n V6 

(T? 2 -l)? 2 

^(v^&i ) when £>2 > 



2 r 1 /-Aqcd 

/l,2(? ) = 4e d 1ta em m D f D / rfx! / ^ 1 fo 1 a 2 (f2)exp[-5 B (?2)]te(xi)- r ^ 

iO Jo {T}- 

f^K^^/Ihbx ) when D 2 > 
^0(^1) when£> 2 <0 



x 2 mD (477 4 - 577 2 - 37] J77 2 - 1 +477 V 77 2 - 1 + 1) +m B (xi - 2)(2t] 3 + 2t} 2 J77 2 - 1 - J77 2 - 1 - 277) , 



r 1 /-Aqcd 

/2.2(? ) = -ledizaemniBfD / dx\ / db\ b\ a 2 (f 2 ) exp[-S B (r 2 )] <j> B (xi ) 

./() .70 (JJ Z 

f i«b(v^i»i) when£» 2 >0 
I ? H o(V\D2~\ b i) when£» 2 <0 



v/6 



(i? 2 -ik 2 

KoiVlhbi ) when Z> 2 > 



x ^ B (x 1 - 2)(77 2 + 77 ^77 2 -l-l) + 2 mD (277 3 +277 2 x /77 2 -l-\/77 2 - I-277) 

2 f 1 /"Aqcd 

/l,3(9 ) = -4e c iza em m D f B / djr 2 / rfZj2^2fl2(f3)exp[-5 D (f3)]fe(x 2 )-^ — — T 

Jo Jo (77 Z -1)<7 Z 

( ^0(^*2) whenD 3 >0 1 
} iH (^\Di\b 2 ) when£> 3 <0 J 



2m B -277 3 - 277 2 v / 77 2 - 1 + Jt7 2 - I+277 



+m D x 2 (477 4 - 577 2 - 3T7 J77 2 - 1 +4T? 3 J77 2 - 1 + 1) - 2(t7 2 + 77 J77 2 -1-1) 



fl /-Aqcd 

f2,3(9) = -4e c 7tOe m m B f B / dx% / ^2^2"2(^)exp[-5z)(f3)]fc(x 2 ) 



X 



"Jo " (?7 2 -l)g 2 
±K (VD^b 2 ) whenD 3 >0 1 
|#o(\/T5#2) when£> 3 < J 



x ^2 mij (77 2 + 77 \/77 2 -l-l)+m D ^ 2 (-277 3 -277 2 \/ 77 2 -l + \/77 2 -l+277) + 2\/77 2 -l 

fl /-Aqcd 1^/3 

fi.4{q ) = -4e u lza em m D f B / dx 2 I db 2 b 2 a 2 {t A )&xp[-S D {t A )]^ D {x 2 )—=— — —^H (^/D 4 b 2 ) 

Jo Jo v2(?7 2 - l)q 2 



x 2 rajJ 27 7 3 +277 2 v / 77 2 -l- J77 2 - 1-2tj +m D (x 2 - 2) U77 4 - 5i7 2 - 3?7 J77 2 - 1 + 417 3 - 1 + 1 



r 1 /-Aqcd ('y^ 

f2MQ ) = -4e«7rOe m m B / B / tfx 2 / ^ 2 i2fl2(U)exp[-5 D (f 4 )]te(x 2 )— =— — ——H (^/D 4 b 2 ) 

Jo Jo V2(f7 2 -l)g 2 



-2m B (77 2 + 77 ^77 2 - 1 - 1) - mo (x 2 - 2) \2rj 3 + 277 2 ^?7 2 - 1 - ,/ 77 2 - 1 - 277 ) ) , (19) 
where (z) = Jq(z) + iY (z), and J , Y and K Q are Bessel functions. 



The expressions forD, (i = 1,2,3,4) are given as 



D\ = -m 2 D + m 2 B + m B m D x\(r\ + \J y] 1 — 1) 



D 2 = -ml(l -xi)-m 2 D - m B m D [-2r] +x\(r] + ^/r\ z - 1)] 



£»3 = -trig + m 2 D + m B m D X2(ri + \/r\ 2 - 1), 



D 4 = -ml - 2m 2 D + m B m D (r] - \J r\ 2 - 1). (20) 

The hard scale f's in the amplitudes are taken as the largest energy scale in the hard kernel Hq (or Ko): tj = 
max ( l/bj^j with j = 1 when i = 1,2 and j = 2 when i = 3,4. 

Functions, S# and result from summing both double logarithms caused by soft gluon corrections and singular 
ones due to the renormalization of ultra-violet divergence. Sb.d are defined as 

S B (t)=s( Xl P+M)+2 f ^r 9 (M')> (2D 

J\/b\ jtl 

S D (t)=s(x 2 P?,b 3 )+2 [' %7 q {li'), (22) 
where s(Q,b), so-called Sudakov factor, is given in |[T0l as 

Q 1 asin') 



Jl/b M 



(23) 



f2, , Q ] a v (A- 

| -(2 7£ - 1 - log2) + C F log ^ | ^ 

f 67 % 2 10 2„, y £ l /a s (jLi')\ 2 , Q 

+ < n f + -Bo oa— ^ og ^ 

\ 9 3 21 1 3 H ° 2 )\ n ) B /x' 

where = 0.57722 • is Euler constant, and y q = a s / % is the quark anomalous dimension. 

References 

[1] M. Antonelli et al, Phys. Rept. 494, 197 (2010) larXiv:0907.5386i rhep-phll. 

[2] W. Altmannshofer, P. Ball, A. Bharucha, A. J. Buras, D. M. Straub and M. Wick, JHEP 0901, 019 (2009) 
HarXiv:081 1.12141 [hep-ph]] . 

[3] A. J. Buras. larXiv: 1 102.56501 rhep-phl . 

[4] M. Beneke, T. Feldmann and D. Seidel, Nucl. Phys. B 612, 25 (2001) ||arXiv:hep-ph/01 06067 1. 



[5] O. Antipin and G. Valencia, Phys. Rev. D 74, 054015 (2006) |arXiv:h ep-ph/0606065|. 



[6] J. A. Macdonald Sorensen and J. O. Eeg, Phys. Rev. D 75, 034015 (2007) [arXiv:hep-ph/0605078|. 

[7] N. Kivel. larXiv:0708.2393l rhep-phl. 

[8] Y. Li and C. D. Lu, Phys. Rev. D 74, 097502 (2006) | |arXiv:hep-ph/0605220| . 

[9] Y. Y. Keum, H. n. Li and A. I. Sanda, Phys. Lett. B 504, 6 (2001) | |arXiv:hep-ph/0004004) . 



[10] Y. Y. Keum, H. N. Li and A. I. Sanda, Phys. Rev. D 63, 054008 (2001) QarXiv:hep -ph/0004 173) ; 
C. D. Lu, K. Ukai and M. Z. Yang, Phys. Rev. D 63, 074009 (2001) | |arXiv:hep-ph/00042 1 3| ; 
A. Ali, G. Kramer, Y. Li, C. D. Lu, Y. L. Shen, W. Wang and Y. M. Wang, Phys. Rev. D 76, 074018 (2007) 
| arXiv:hep-ph/0703162| . 

[11] S. Descotes-Genon and C. T. Sachrajda, Nucl. Phys. B 625, 239 (2002) larXiv:hep-ph/0109260| . 

[12] F. Feng, J. P. Ma and Q. Wang, Phys. Lett. B 674, 176 (2009) llarXiv:0807.0296i rhep-phll; 
H. n. Li and S. Mishima, Phys. Lett. B 674, 182 (2009) HarXiv:0808. 1 5261 [hep-ph] 1 ; 
F. Feng, J. P. Ma and Q. Wang, Phys. Lett. B 677, 121 (2009) narXiv:0808.4017i rhep-phll. 

[13] H. n. Li and S. Mishima, Phys. Rev. D 71, 054025 (2005) larXiv:hep-ph/0411146| ; 

H. n. Li, S. Mishima and A. I. Sanda, Phys. Rev. D 72, 1 14005 (2005) | |arXiv:hep-ph/05080411 ; 
H. n. Li and S. Mishima, Phys. Rev. D 74, 094020 (2006) ( arXiv:hep-ph/0608277| . 

[14] G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. Mod. Phys. 68, 1 125 (1996) flarXiv:hep-ph/95 123801 . 

[15] R. H. Li, C. D. Lu and H. Zou, Phys. Rev. D 78, 014018 (2008) HarXiv:0803.1073ll [hep-phll. 



9 



